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Linearization of coupled second order nonlinear ordinary differential equations (SNODEs) 
is one of the open and challenging problems in the theory of differential equations. 
In this paper we describe a simple and straightforward method to derive linearizing 
transformations for a class of two coupled SNODEs. Our procedure gives several 
new types of linearizing transformations of both invertible and non-invertible kinds. 
In both the cases we provide algorithms to derive the general solution of the given 
SNODE. We illustrate the theory with potentially important examples. 
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1. Introduction 



Continuing our study on the integrability and linearization of coupled second order 
nonlinear ordinary differential equations (SNODEs), in the present paper, we focus 
our attention on the linearization of two coupled SNODEs. The present study arises 
not only for the completeness of Part IV (Chandrasekar et al. 2009) but also to show 
the importance of unfinished tasks that exist in the theory of linearization of two 
coupled SNODEs. As far as the first point is concerned, we show that one can also 
solve a class of coupled second order nonlinear ODEs by transforming them into 
two second order free particle equations and from the solutions of the latter one 
can construct the solution of the former, even though this is a nontrivial problem 
in many situations (one can also transform coupled nonlinear ODEs into uncoupled 
nonlinear ones, which has already been pointed out by us in the previous paper, 
Part IV). Regarding the second point, we wish to stress the fact that linearization 
of coupled nonlinear ODEs is a vast area of research which is still in its early 
stage. In this paper, we show that in spite of the difficulties that exist in this topic, 
one can make useful progress on certain issues, namely (i) developing a method 
to deduce all linearizing transformations wherein the new dependent variables are 
functions of the old dependent and independent variables only and not derivatives 
of the dependent variables, and (ii) developing a method of constructing solutions 
of nonlinear ODEs from the linear ones in the case of non-point transformations. 
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Even though the modern theory of linearization of nonlinear ODEs had origi- 
nated and developed with the works of Lie, Tresse and Cartan (Mahomed & Leach 
1989; Steeb 1993; Olver 1995; Ibragimov 1999; Chandrasekar et al. 2005), the en- 
tire subject was lying dormant for more than a century. Only recently, during 
the past two decades or so, notable progress has been made to linearize nonlin- 
ear ODEs through non-point (Duarte et al. 1994) or generalized transformations 
(Chandrasekar et al. 2006). For example, focussing our attention on single second 
order ODEs, generalized Sundmann (Eulcr et al. 2003; Eulcr et al. 2004) and gener- 
alized linearizing transformations (Chandrasekar et al. 2006) have been introduced 
to linearize a class of equations that cannot be linearized by invertible point trans- 
formations. As far as two coupled SNODEs are concerned, to our knowledge, most 
of the studies were focussed only on invertible point transformations irrespective of 
whether it is an analytical approach or a geometrical formulation. For a survey on 
this topic one may refer the recent papers of Merker (2006) and Mahomed & Qadir 
(2007) and for the earlier works in this direction we cite Crampin et al. (1996), Fels 
(1995), Grossman (2000), Soh & Mahomed (2001) and Qadir (2007). More recently 
Sookmee & Meleshko (2008) proposed a new algorithm to linearize the coupled 
second order ODEs by sequentially reducing the order of the equation. 

In this work we aim to give a new dimension to the theoretical development 
of linearization of nonlinear dynamical systems having two degrees of freedom by 
proving that one can unearth a wide class of linearizing transformations besides 
invertible point transformations. Of course, the latter ones form a sub-class of the 
new ones which we construct in this paper. In the present study, we not only derive 
several new types of linearizing transformations but also propose systematic proce- 
dures to derive the general solution in all these cases. We also wish to emphasize 
here that we derive all these transformations from the first two integrals alone and 
thereby establish a potentially simple, straightforward and powerful approach in 
the theory of differential equations. 

The plan of the paper is as follows. In §2 a we briefly describe the method of 
deriving linearizing transformations for a system of two coupled second-order ODEs. 
We show that one can have two classes of linearizing transformations, depending 
upon the nature of the independent variables. If the new independent variables are 
the same (zi — Z2), we put them in Class-A category, and if they are different (z\ 7^ 
22) j then we put them in Class-B category. In §2 b(i) we consider Class-A category 
and identify three types of linearizing transformations. In §2 b(ii) we consider Class- 
B category and identify six types of linearizing transformations. In §3 we consider 
one specific example for each of the nine types of linearizing transformations we 
have identified, and obtain general solution to each one of them to demonstrate our 
procedure. Finally, we present our conclusions in §4. 

2. Linearizing transformations 

(a) Method of deriving linearizing transformations 

To begin with let us consider a system of two coupled second order nonlinear 
ODEs, R{t,x}, (vide equation (2.1) in Part IV (Chandrasekar et al. 2009)), 

x = <t>i(t,x,y,x,y), and y = 4> 2 {t, x,y, x,y). (2.1) 
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Any transformation of the form T{t,x} defined by 

wi = fi{t,x,y), zi = J f 3 (t,x,y,x,y)dt, 

w 2 = f2{t,x,y), z 2 = J fi(t,x,y,x,y)dt, (2.2) 

which transforms the given set of nonlinear ODEs (2.1) to the free particle equa- 
tions, 



d 2 wx _ d 2 w 2 
dz\ dz 2 



2 -o. -^r = ' ( 2 -3) 



is called a linearizing transformation in the present work. 
Let 

h = T(t,x,y,x,y), and I 2 = Q(t,x,y,x,y). (2.4) 

be the two first integrals of motion of the coupled system (2.1) and that they can 
be explicitly found, if they exist, for example, by using generalized modified Prelle- 
Singcr (PS) method formulated in Paper IV (Chandrasckar et al. 2009). Then the 
following theorem ensures that the transformation can be deduced from the /j's, 
i = 1,2. 

Theorem: Suppose a given nonlinear system R{t,x} of ODEs (2.1) is linearizable 
to a system of two uncoupled free particle equations through the linearizing trans- 
formation T{t,x} of the form (2.2), then the latter can be deduced from the first 
integrals Ij(t, x, y, x, y), i = 1,2. 

Proof: Let us re-express each of the functions T and Q in (2.4) as a product of two 
new functions, that is, 

h = -tt: — 1 . 7z ^rfi(t,x,y), 
h(t,x,y,x,y) dt 

12 h(t,x, y,x,y) dt^ t,X,V ^' ^'^ 

Again rewriting / 3 and f± as total time derivatives of another set of functions, say 
z\ and z 2 , respectively, that is dz\/dt = fo{t, x, y, x, y) and dz 2 /dt = fi{t, x, y, x, y), 
equation (2.5) can be further recast as 



1 dh dh T 1 df 2 df 



i-fydt-dzS 2 -^dt-dz 2 - 

Now identifying the functions fi(t,x,y) = w\ and f 2 (t,x,y) — w 2 as the new 
dependent variables, equation (2.6) can be further recast in the form 

<*u>i = f dw 2 _ f , 2 J , 

dzi ' dz 2 

where 1\ and I 2 are the redefined constants. Obviously, (2.3) follows straightfor- 
wardly from (2.7). Consequently the new variables, Zi and u>i, i = 1,2, defined by 
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equation (2.2) help us to transform the given set of coupled SNODEs into two linear 
second order ODEs which in turn lead to the solution by trivial integration. The 
variables w^s and z^s, i = 1,2, then define the linearizing transformations for the 
given equation (2.1). □ 
It may be noted that in general the new dependent variables, Wi and w 2 , may 
also involve x and y, that is w\ = Gi(t,x,y,x,y) and w 2 = G 3 (t,x,y,x,y), and 
this possibility may lead us to identify more generalized transformations like point- 
contact and generalized-contact transformations. However, in this paper we will 
confine ourselves only to the forms of w\ and w 2 given by (2.2). 

(b) The nature of transformations 

An important question which we will focus upon in the present paper is that 
what are the possible forms of linearizing transformations one can unearth through 
the above procedure. We recall here that in the case of scalar SNODEs, one has 
point, generalized Sundman and generalized linearizing transformations (Chan- 
drasekar et al. 2005;2006). As far as the coupled SNODEs (2.1) are concerned, 
since there are two independent variables Z\ and z 2 as given in equation (2.2), one 
can choose them to be either the same, Z\ = z 2 (Class-A), or different, Z\ ^ z 2 
(Class-B). In the case of Class-A transformations, one can construct three different 
types of linearizing transformations, while for class-B one can formulate six differ- 
ent types of linearizing transformations, as we point out below. However, we also 
note that even further types of local transformations involving the variables x and 
y are possible but these are not included in the present study. 

(i) Class-A Linearizing Transformations (z\ = z 2 = z) 

In the case of Class A transformations we have w\ = fi(t, x, y), w 2 — f 2 (t, x, y), 
z\ = z 2 = z = J f 3 (t, x, y, x, y)dt = J f±(t, x, y, x, y)dt. Now appropriately restrict- 
ing the form of fa (= one can identify three different types of linearizing 
transformations. 

1. Suppose z\ = z 2 = z is a perfect differential function and w^s, i = 1,2, and 
z do not contain the variables x and y, then we call the resultant transfor- 
mation, namely, Wi = fi{t,x,y), w 2 = f 2 (t,x,y), z = f 3 (t,x,y), as a point 
transformation of type-I. 

2. On the other hand, if z is not a perfect differential function, and w^'s, i = 
1,2, and z do not contain the variables x and y, then we call the resultant 
transformation, namely, wi = fi(t,x,y), w 2 = f 2 (t,x,y), z = J f 3 (t 7 x,y)dt 7 
as a generalized Sundman transformation of type-I. 

3. As a more general case, if we consider the independent variable z to con- 
tain the derivative terms also, that is wi = f\(t,x,y), w 2 — f 2 (t,x,y), 
z = J f 3 (t,x,y,x,y)dt, then we call the resultant transformation as a gen- 
eralized linearizing transformation of type-I. 

In our analysis, we do not consider the possibility ui\ = fi(t,x,y), w 2 = f 2 (t,x,y) 
and z = f 3 (t 7 x,y,x,y) because the procedure to handle it is different from the 
presently discussed linearizing transformations. This possibility will be studied sep- 
arately. 
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(ii) Class- B Linearizing Transformations (z\ 7^ z 2 ) 

In the Class-B type of linearizing transformations we have w\ — f\(t,x,y), w 2 = 
f 2 (t,x,y) and z 1 = J f 3 (t,x,y,x,y)dt, z 2 = J f 4 (t,x,y,x,y)dt, z 1 ^ z 2 . Now ap- 
propriately restricting the forms of f 3 and f 4 , one can get six different types of 
linearizing transformations. 

1. If zi and z 2 are perfect differential functions and w^s and Zi's, i = 1,2, do 
not contain the variables x and y, then we call the resultant transformation, 
namely, ivx = h(t,x,y), w 2 = f 2 (t,x,y), z x = f 3 (t,x,y), z 2 = f 4 (t,x,y), as a 
point transformation of type-II. 

2. Suppose Z\ is a perfect differential function and z 2 is not a perfect differential 
function or vice versa, and if z\ and z 2 do not contain the variables x and 
y, then we can call the resultant transformation, namely wi = fi(t,x,y), 
w 2 = f 2 (t,x,y), z\ = f 3 (t,x,y), z 2 = J f 4 (t,x,y)dt or z x = J f 3 (t,x,y)dt, 
z 2 = f 4 (t, x, y), as a mixed point-generalized Sundman transformation. 

3. On the other hand, if any one of the independent variables contains the 
variables x and y, then we call the resultant transformation, namely wi — 
fi(t,x,y), w 2 = f 2 (t,x,y), z\ = f 3 (t,x,y), z 2 = J f 4 (t,x,y,x,y)dt or z\ = 
/ f?>{iTX,y,x,y)dt and z 2 = f4(t,x,y), as a mixed point-generalized lineariz- 
ing transformation. 

4. Suppose the independent variables are not perfect differential functions and 
are also not functions of x and y, that is, w\ = fi(t,x,y),w 2 = f 2 {t,x,y), z\ — 
J f 3 (t,x,y)dt, z 2 = J fi{t,x,y)dt, then we call the resultant transformation 
as a generalized Sundman transformation of type-II. 

5. Further, if one of the independent variables, say z\, does not contain the 
derivative terms while the other independent variable z 2 does contain the 
derivative terms or vice versa, that is w\ = fi(t,x,y), w 2 = f 2 (t,x,y), z\ = 
I h(t,x,y)dt, z 2 = J f 4 (t,x,y,x,y)dt or z\ = J f 3 (t,x,y,x,y)dt, z 2 = J f 4 (t, 
x,y)dt, then we call the resultant transformation as a mixed generalized 
Sundman-generalized linearizing transformation. 

6. As a general case, if we allow both the independent variables, z\ and z 2l 
to be non-perfect differential functions and also to contain derivative terms, 
that is, wi = h(t,x,y), w 2 = f 2 (t,x,y), zi = J f 3 (t,x,y,x,y)dt, z 2 = 
J f 4 (t, x, y, x, y)dt, then the resultant transformation will be termed as a gen- 
eralized linearizing transformation of type-II. 

Finally the possibility that wi = w 2 = f 2 (t,x,y), z\ = f 3 (t,x,y,x,y) 

and z 2 = f 4 (t 7 x, y, x, y) is not considered in the present study and will be pursued 
separately. 

3. Applications 

In this section, we consider specific examples and illustrate each one of the lin- 
earizing transformations identified in the previous section so as to make clear the 
applicability of them under different situations. 
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(a) Class- A Linearizing Transformations (z\ = z 2 ) 
(i) Example: 1 Point transformation of Type-I 
Let us consider the system of SNODEs, 

* + fc^+^ = 0, g-fc^+^o, (3.1) 
2xy(x - y) 2xy(x - y) 

where u is a arbitrary constant. The two first integrals associated with (3.1), which 
can be obtained using the formulation given in §2 of Part IV (Chandrasekar et al. 
2009), can be written as 

1\ = (x + y) sin(wt) — lo(x + y) cos(wt), 

I 2 = ((xy + xy) sin(ii>£) — 2ujxy cos(bjt)). (3.2) 



Rewriting (3.2) in the form (2.5)-(2.6), we get 



. 2 . d .. . . ,.. dt dwi dwi 

h = sin (ut) — ((x + y)cosec{ut)) = — = — - , 

at az\ dt dz\ 

.2/ n d / , , u dt dwi dw 2 ,„ „. 

h = sin M^(^ecM)) = — — = — , (3.3) 

so that we can identify point transformation of type I as 

w\ = (x + ?y)coscc(o)t), w 2 = y^coseciujt) , z\ = z 2 = z = — — cot(wt). (3.4) 

U) 

Utilizing the transformation (3.4), one can transform (3.1) to a set of free particle 
.2 ,2 

equations, namely — and -^f- = 0, so that w\ = I\z + I^ and w 2 = I 2 z + l4, 
where ^3 and I4 are integration constants. Substituting the expressions for w^'s, 
i = 1,2, and z in the free particle solutions and rewriting the resultant expressions 
in terms of the old variables x and y, one obtains the general solution for equation 
(3.1) in the form 

x(t) = -±(A T VA2-4B2), y(t) = ±(A±y/A2-4B*), (3.5) 

where A = I\ cos(wt) + 7 3 sin(a>£) and B — I 2 cos(wt) + I4 sin(wi). Here we point 
out that the nonlinear system (3.1) admits amplitude independent frequency of 
oscillations. 

In the above example we have considered the new dependent variables (w\, w 2 ) 
and independent variable z to be functions of x, y and t only. We will now consider 
examples which admit more general transformations. 

(ii) Example:2 Generalized Sundman Transformation of Type-I 

Let us focus our attention on the two-dimensional Mathews-Lakshmanan oscil- 
lator system of the form (Carinena et al. 2004; Chandrasekar et al. 2009) 

\(x 2 + y 2 + X(xy — yx) 2 )x — a 2 x 



x = 



(1 + Ar 2 ) 



.. _ A(i: 2 + y 2 + \{xy - yx) 2 )y - a 2 y 

y- (i + xr 2 ) ' (6 - b) 
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where r — \J x 1 + y 2 and A and a are arbitrary parameters. For a = 0, equation 
(3.6) admits the following two integrals of motion, 

f _ (1 + Ay 2 )± - Xxyy f _ (1 + Aa; 2 )y - Axyi 

We note that the integrals I 2 and I3 (vide equations (5.35) and (5.36) in Part 
IV (Chandrasekar et al. 2009)) can be derived from (3.7) by using the relations 
h = + If + ^i) and I 3 = (h + I 2 ) 2 - 2A/ 2 , where h is given in (5.35) 

of Part IV. The general case (a / 0), can be linearized through mixed generalized 
Sundman - generalized linearized transformation (see example 8 below) . To demon- 
strate the linearization through generalized Sundman transformation of Type-I, we 
here consider equation (3.6) with a = 0. 

For the case a = 0, one may note that on making a substitution y(t) = y(x{t)) 
into (3.6), one can obtain a non-autonomous second order ODE in y(x). Though 
this equation satisfies the linearization condition for point transformation (Sookmec 
& Meleshko 2008), finding the linearizing transformation and the general solution 
for the transformed ODE turns out to be nontrivial. On the other hand, we provide 
a straightforward procedure of linearization. 

The above two integrals (3.7) can be rewritten as 

t , , 9n d ( x \ dt dwi dwi 

h ^ 1 + A ^(7TTP)^^ = V (3 ' 8) 

e^VVl + Ar 2 / dz 2 dt dz 2 V ' 

From the above equation, we identify the new dependent and independent variables 
as 



y f dt 

«* = ^™*. W2 = 7tt^' Z1 = Z2 = Z = J JTT^-y (3 - 10) 



X 

VlTAr 2 : 



One may observe that the independent variables z\ and z 2 are not perfect differen- 
tials even though they turn out to be identical. By utilizing the above new variables, 
one can transform equation (3.6), with a = 0, to the free particle equations, that 
is,^ = 0and^ = 0. 

Unlike the earlier example, one cannot unambiguously integrate these two lin- 
ear equations in terms of the original variables due to the nonlocal nature of the 
independent variable. To overcome this problem one should express (1 + \r 2 ) in 
terms of either z\ or z 2 so that the resultant expression dz\ = dz 2 = ^x r ' 2 ) can 
be integrated. In the following, we describe a procedure to obtain an expression for 
the new independent variable. 

Now integrating equation (3.8) and using first relation in (3.10), we get 

hzi, (3.11) 



Vl + Ar 2 

where we have fixed the integration constant to be zero (without loss of generality). 

dwi I± 

dw 2 f 2 



On the other hand from expressions (3.8) and (3.9) we get = y- from which 
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one obtains 



h 



w 2 + h 



h 



y 



h ' - VI + Ar 2 I 2 VI + Ar 2 

where I3 is the integration constant. Equation (3.12) provides us an identity 



(3.12) 



VI + Ar 2 " h 
Now squaring and adding the equations (3.11) and (3.13), we obtain 



\ r 2 f 2 f f 2 f 2 



A(W + ^K 



1 + Ar 2 ^ ' / : ' If 

From (3.14), one can express (1 + Ar 2 ) in terms of z\ as 

1 + Ar 2 = 



)■ 



1 



l-A((7 1 2 +/I), 1 2 -^i + 7?7? 



/1 



(3.13) 



(3.14) 



(3.15) 



Substituting (3.15) in the last relation given in equation (3.10), we arrive at the 
following integral relationship between z\ and t, namely, 



dzi 



i-A( (I't + iDzf 



'2 J 3 



r2 r2 



(3.16) 



'1 " 7 i 

Since the variables are separated now, one can integrate this equation and obtain 
an expression which connects the new independent variable with the old indepen- 
dent variable through the relation 



Z\ 



VA/2/3 - iiujta,n[uj(t - to)] 



(/ 2 (A/ 2 - 1) - 7 2 ) > 0, 



(3.17) 



Va/i(/ 1 2 + 4 2 ) 

where uj = V\\J Jf (A/| — 1) — If and t is the fourth integration constant. From 
equations (3.11)-(3.13) and (3.15), we get 



x(t) = I\Z\ 



1-X (If+Ii)z(- 



'2 J 3 



h 



-Zl 



y(t) = hzx 



h 



l-\Uli + If)zi- 



2IU. 



h 



-zi 



il 



.(3.18) 



Substituting the expression (3.17) in the above equation (3.18) and simplifying 
the resultant expressions, we arrive at the following general solution for equation 
(3.6) with a = in the form 



(3.19) 



x(t) = A yXllh cos[u(i - t )} - u>h sm[w(t - to)] ) , 
y(t) = — IzAi XI1I3 cos[w(i — to)] + ujsin[Lu(t — t 



where A = _i_./_JS±iaL an d w = vW/f (A/ 2 - 1) - If. 
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(iii) Example:3 Generalized linearizing transformation of Type-I 

In the previous example we restricted the new independent variable to be a 
nonlocal one and be a function of t, x and y only. Now we relax the latter condition 
and also allow the independent variables Z\ and z 2 to contain derivative terms, 
namely x and y. To illustrate this case let us consider a two coupled second order 
equations of the form 



x = 



2xy(xx + x 3 ) — 2yx 3 
x^y 

2xy(x 2 y — yx + xy) + x(x 2 y 2 — y 2 x 2 ) 

y 3 * 

x°y 

One can easily identify two integrals for the equation (3.20) in the form 



h 



1 1 

y\ — + - 



and 



, y v 
h = y - - -r 



(3.20) 



(3.21) 



Rewriting the above integrals as 



y 



^ x dt 



1 



dwi 
dz\ ' 



h = — — 

x dt 



(log 


X 






) 




.y. 





= (3.22) 

dz 2 



we identify the following set of linearizing transformations for the equation (3.20), 
that is 



Wl = t — 



w 2 = log 



Z\ = z 2 



rdt. 



y 



(3.23) 



One may note that the independent variables are not only nonlocal but also involve 
derivative terms. It is easy to check that equation (3.23) transforms (3.20) to the 
linearized form (2.3). 

Again, as in the previous example, one cannot directly obtain the solution for 
(3.20) by direct integration of the linear ODEs due to the nonlocal nature of the 
independent variables. This can be overcome by expressing the term x/y 2 in terms 
of either z\ or z 2 so that the resultant equation can be integrated to obtain an 
explicit form for the new independent variable in terms of the old variables as we 
discuss below. 

From (3.21), we have 



x 

7/2 



1 



(3.24) 



Since = I\ (from equation (3.22)), we have wi — 1\Z\, so that 



t = hz-L, 

x 



(3.25) 



where without loss of generality we have fixed the integration constant to be zero. 

h 



On the other hand, from (3.22), we have — j^, which in turn gives 



J l ^ T 
^2 



(3.26) 
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In terms of old variables (3.26) reads as 

+ h- (3.27) 



* 1 h i 
t = — log 

x h 



From the identites (3.25) and (3.27), we can express y/x interms of z\ in the form 

|=exp[-^(Ji«i-/ 3 )]. (3.28) 

Now substituting (3.28) into equation (3.24), we can express (p-) in terms of Z\, 
and plugging the latter relation into the third relation in (3.23), we arrive at 

dt = (h - exp[-2^(/ 1 z 1 - h)])d Zl . (3.29) 

Integrating the above equation, we obtain 

e -^(i lZl -i 3 ) 

t + to = hzi + - 7T T , (3-30) 

t— - 

where to is the fourth integration constant. Substituting the expression z\ = —pf- 
(vide equation (3.25)) into equations (3.28) and (3.30), we obtain the general solu- 
tion of (3.20) in the implicit form 

a! ( t ) = -- + -y___ , y(t) = x(t)eM-f i (t- W) -Is)]. (3.31) 

(b) Class-B Linearizing Transformations (z\ ^ z 2 ) 

In the Class-A category, in all the three examples, we considered that the new 
independent variables, z\ and Z2, are identical. However, this need not always be 
the case in the theory of linearizing transfomations, as discussed in Sec. §2. We now 
present specific examples to illustrate more general transformations. 

(i) Example:^ Point transformation of Type-II 

Let us consider a quasi periodic oscillator governed by a set of two coupled 
SNODEs of the form 

(xy- yx) 2 + 2x 2 y(ujl(x + y) - 2u\y) = Q 
2xy(x-y) 

~ _ {xy - yx) 2 + 2xy 2 (uj 2 (x + y) - 2u%x) =Q 
2xy(x — y) 

To explore the linearizing transformation for equation (3.32), we consider the 
two associated integrals, 

l\ = (x + y) coswit + u>i(x + y) smu)-\t, 

I2 = — — {{xy + yx) smui2t — 2u>2xy coso^i)- (3.33) 
2Jxy 
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j2 



Rewriting (3.33) in the form 

h = cos 2 {Lo 1 t)j t {{x + y)scc{Lo 1 t)) = (3.34) 

7 2 = s\n 2 {Lo 2 t)^{^/xycosec{uj2t)) = ^j^-, (3.35) 
one can identify the new dependent and the independent variables as 
wi = (x + y) scc(wii), z\ = — tan(u>ii), 

LUX 

w 2 = \fxy 'COsec(a;2^), %2 = cot(uj2t). (3.36) 

UJl 

One may note that now the independent variables z\ and Z2 are not the same. 

,2 

The new variables transform (3.32) to the free particle equations -^r = and 

0. From the general solutions w\ = I\Z\ + I3 and u>2 = I2Z2 + I4, where 
Ii's, i — 1,2,3,4, are the integration constants, and using the expressions for Wi's 
and ZiS, i — 1,2, given in (3.36), we arrive at the general solution for the equation 
(3.32) in the form 

x{t) = ~(A T VA 2 -4B*), y(t) = ±(A±VA 2 -4B 2 ), (3.37) 
where now A~l\ sin(wit) + J3 cos(wit) and B = h cos(uj2t) + h sm(uj2t)- 

(ii) Example:5 Mixed Point-Generalized Sundman transformation (PGST) 

Let us consider the two dimensional force-free coupled Duffing- van der Pol os- 
cillator (DVP) equation of the form 

x + 4(a + (3(kix + k 2 y) 2 )x + a(3a + A[3{k x x + k 2 y) 2 )x = 0, 

y + 4(a + (3{k lX + k 2 y) 2 )y + a(3a + A(3{k lX + k 2 y) 2 )y = 0. (3.38) 

One may note that the point transformation X = k\x + k 2 y and Y = k\x — k 2 y 
helps one to rewrite (3.38) in a separable form, 

X + 4(a + f)X 2 )X + a(3a + Af3X 2 )X = 0, (3.39a) 
Y + 4(a + (3X 2 )Y + a(3a + 4(3X 2 )Y = 0. (3.39b) 

The solution to equation (3.39a) can only be obtained in implicit form (Chan- 
drasekar et al. 2005). Consequently, equation (3.39b) cannot be solved explicitly in 
this way. Further, the linearization of the scalar DVP oscillator (3.39a) itself has 
not yet been reported. In the following we use our procedure to find the linearizing 
transformation and general solution to (3.38) straightforwardly. 

The first two integrals for (3.38) can be easily identified using the procedure 
given in Chandrasekar et al. (2009) in the form 



y + ay 



h = (k 1 x + k 2 y + a(k 1 x + k2y) + ^-(k 1 x + k 2 y) 3 )e 3at . (3.40) 
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The above integrals can be rewritten as 

y 2 e at d x dw\ 
h = ~ 77 : : n ~(~) 



(y + ay) dt y dz\ 



(3.41) 



-{{k lX + k 2 y)^ at f 



dt 



-2at 



dwi 
dz 2 ' 



from which we can obtain the following linearizing transformations, 

x 



Wi = 



Z\ 



y 



W2 = (l(k 1 x + k 2 y)- 2 + ^)e-^\ 
2 6a 



-at 



v 



Z2 



I 



((k lX + k 2 y)e^ at )- 3 dt- 



(3.42) 



One may note that in the present problem one of the new independent variables, 
that is, z 2 is in a nonlocal form. In terms of the above new variables, equation (3.38) 
assumes the forms ^ = and = 0. 

dzf dz^ 

Now we seek the general solution of (3.38) from the linearized equations. Inte- 

.2 

grating = 0, we obtain 



toi = hzx + h 



(3.43) 



where I3 is the integration constant. Rewriting (3.43) in terms of the old variables, 
we get 



x = he- at + I 3 y. 



(3.44) 



,2 

However, the second linear equation, -^f- — 0, cannot be integrated straightfor- 
wardly (in terms of the original variables) due to the nonlocal nature of the second 
independent variable. To obtain an explicit form for z 2 , we rewrite I 2 in the integral 
form (vide equation (3.41)) to obtain 

{\(k lX + k 2 y)- 2 + ^)e- 2at = -h J [(hx + k 2 y)el at ]- 3 dt = I 2 z 2 . (3.45) 

Equation (3.45) provides us a relation 



{hx + hy)- 1 = \j 2I 2 z 2 e 2at 



M 
3a' 



(3.46) 



Now substituting the relation (3.46) into the non-local variable z 2 (vide equation 
(3.42)), one gets 



— - -I 2 {2I 2 z 2 - - 
Solving the above equation, we obtain 



t -^e 2at = — 



= -h{2hz 2 - ^-e- 2at )^e~ 2at . 

V3 



2a 



3/2 



V3 tan 



\2a{2I 2 z 2 -^e- 2 ^)h -1 
(l + a(2I 2 z 2 -^e- 2at )?) 2 



(3.47) 



(3.48) 



1 - a(2I 2 z 2 - ^-e~ 2at )h + a 2 {2I 2 z 2 e 2at - ^- e - 2at ) 
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where a = y — irj and to is the fourth integration constant. From the expression 
(3.48) and equations (3.44) and (3.46) one can deduce the general solution for 
(3.38) in implicit form. The resultant expression coincides exactly with equation 
(6.18) given in Chandrasekar et al. (2009). 

In the present example we considered one of the independent variables to be 
in a nonlocal form. As we have two independent variables one can also have the 
possibility of having both the independent variables to be of nonlocal form. Indeed 
this is the case of our next example. 



(iii) Example:6 Generalized Sundman transformation of Type-II ( GST-II) 
To illustrate the GST-II, we consider the equation of the form 
2 2 

' x - ( x 2 + y 2^ ±2 - y 2 ) x + 2 ^y) + p x = °' 

The first two integrals for (3.49) can be evaluated as 

(2(ty - y)xy) + (x 2 - y 2 ){tx - x) 



h 



t 2 (x 2 + y 2 ) 2 



T _ (2xy(tx - x) + {x 2 - y 2 )(ty - y)) 

h - t 2 (x 2 + y 2 ) 2 • (3 - 50) 



Rewriting these two integrals as 



x 2 d { x y 2 \ dwi 

y 2 d ( y x 2 \ dw 2 , 



{x 2 + y 2 ) 2 dt\t ty J dz 2 
we identify the linearizing transformations in a more general form 



x y 2 y x 2 ({x 2 + y 2 ) 2 rU f(x 2 + y 2 ) 2 
l + Tx> W2 = -t + ^ z ^J—^ dt > Z2 ^J^^ ( 



The GST-II, (3.52), takes equation (3.38) to the free particle equations, = 

,2 1 

and -^r — 0. To obtain the solution in terms of the original variables, we have to 

replace both J ((x 2 + y 2 ) 2 /x 2 )dt and / {{x 2 + y 2 ) 2 /y 2 )dt by the variables z\ and t, 
and z 2 and t, respectively, and integrate the resultant equations. 

To do so, first we rewrite the first integrals I\ and I 2 given by (3.51) in integral 
forms, which in turn lead us to w\ — I\Z\ and w 2 = I 2 Z 2- Since W\ and w 2 do not 
contain nonlocal variables we can replace them by the old variables (vide equation 
(3.52)), that is 
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where we have fixed the integration constants to be zero (without loss of generality). 

We observe that to integrate the last two expressions in (3.52) one should further 
replace Z\ and z 2 in terms of t. So we substitute the above expressions for x and y 
in terms of Z\ and z 2 , respectively, in the last two relations in (3.52), and obtain 

dz x = l\z\t 2 dt and dz 2 = llz 2 2 t 2 dt. (3.54) 

Now integrating both the equations, we get 

Zl = (J?(/ 3 - I))" 1 and Z2 = (ll( h -t))-\ (3.55) 

where ^3 and J 4 are the third and forth integration constants, respectively. Plug- 
ging equation (3.55) into (3.53), we arrive at the following general solution for the 
equation (3.49), 



( ) = 3f(Ji(3/ 3 - (If - Ij)t 3 ) - J 2 (3J 4 + 2hl 2 t 

(3/3 - (il - ii)t 3 ) 2 + (3/ 4 + 2W 3 ) 



2 



(t) = 3t(/ 2 ((J 1 2 - - 3/3) - Ji(3J 4 + 2hl 2 t 3 ) 
(3/ 3 - (A 2 - I 2 2 )t 3 ) 2 + (3/ 4 + 2hht 3 ) 2 ' 

where 7 3 = (/1 J 3 - hh)/(ll + h%) and J 4 = (7 2 7 3 + /i/ 4 )/(/ 2 3 - l\h)- 

In the previous two examples we focussed our attention on the case in which 
the new independent variable(s) is(are) nonlocal and does(do) not admit velocity 
dependent terms. Now we relax this condition and allow either one or both the 
independent variables to admit velocity dependent terms but in nonlocal form. 

(iv) Example:! Mixed Point- Generalized linearizing transformation 

To demonstrate this, we consider a variant of the two-dimensional Mathews and 
Lakshmanan equation (3.6) of the form 

.. _ X(x 2 +y 2 + 2X(y-x) 2 )-a 2 .. _ X(± 2 + y 2 + 2X(y - xf ) - a 2 

(l + 2X(x + y)) ' V (l + 2X(x + y)) ' [ ^'> 

Equation (3.57) admits the following two integrals of motion, 

»-*-*> ^-T+^r 3 *' 

Rewriting (3.58) in the form 

h = J t (x-V), (3-59) 



a 2 - A((l + 2X)(y - x) 2 + 2xy) d 
h = 



log[l + 2X{x + y)\ 



(3.60) 



2X{x + y) dt 
one can easily identify the linearizing transformations for (3.57) as 
wi = (x-y), w 2 = log[l + 2A(x + ?/)], 

Zl=t, Z 2 = [ T77 2X }\,t-^ - no TTTTrdt. (3.61) 

J a 2 - X((l + 2X){y -x) 2 + 2xy) V ; 
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In terms of the above new variables, equation (3.57) gets transformed to the free 
particle equations (2.3). One may note that one of the new independent variables 
is not only in non-local form but also contains derivative terms which in turn 
complicates the situation to obtain the general solution. 

Since W\ and Z\ are both of point transformation types, one can integrate the 

.2 

first free particle equation, namely, -^r = and obtain w\ = I\Z\ +I3, where J 3 is 
an integration constant. Replacing the latter in terms of the old variables, one gets 
the relation (x — y) = I\t + I 3 . On the other hand, from the solution of the second 
linear equation = 0, we can write w 2 = I 2 z 2 => log[l + 2X(x + y)] = I 2 z 2 

(again we assume the integration constant to be zero without loss of generality). 

To evaluate z 2 , let us first substitute (3.58) into (3.61), and rewrite the latter 
in the form 

2X{x + y) 2X(h + 2y) 

dZ > = I 2 (l + 2X( X + y)) dt = I 2 (l + 2X( X + y) ) dt - (3 ' 62) 

Now subsituting the form of y (vide equation (3.58)), that is, 

V = ^ ( - A/i ± v /2Aa 2 - A 2 (l + 4A)/? - 2A/ 2 (1 + 2X(x + y))^j , (3.63) 
into equation (3.62), and using the relation (1 + 2X(x + yj) = e l2Z2 , we obtain 

dz 2 = WM-Kl + Wi-U^a (3 . 64) 
Integrating (3.64), we get 

, 2 = i- loe f 2 ° 2 -^+^-"^('-'°) 2 Y ,3.65) 
h \ 2I 2 J 

where to is an integration constant. Substituting the expression (3.65) into the 
relation 2A(x + y) = e l2Z2 — 1, we get 

x + y ^-^ + mi-JilKt-t ?-2 h _ (3 66) 

From equation (3.66), and the relation (x — y) = I\t + I3, we obtain the general 
solution for (3.57) in the form 

_ 2c? - A(l + 4A)Z 2 - 4/|A(t - t ) 2 - Z 2 (2 - AXjht + J 3 )) 
X[ ) ~ 8XI 2 

y{t) = 2a 2 - A(l + 4A)/ 2 - 4/|A(t - t ) 2 - Z 2 (2 + AXjht + J 3 )) _ (3 6?) 

8A/2 

In this example we considered the case in which only one of the independent vari- 
ables is in nonlocal form. Now we consider the case in which both the independent 
variables are in nonlocal forms. 
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(v) Example:8 Mixed Generalized Sundman-Generalized linearizing transformation 

To illustrate this type of linearizing transformation let us consider again equa- 
tion (3.6), but now with a^O. Equation (3.6) admits the following two integrals 
of motion, 



h = {yx - xy), 



(a 2 - X(x 2 + y 2 + X(yx - xy) 2 )) 



1 + Xr 2 



Rewriting these two integrals in the form 



h = y 2 



dt\y 



dwi 



h = 



{a 2 - X{x 2 + y 2 + X(yx - xy) 2 ) d 



2X(xx + yy) dt 
and identifying the linearizing transformations, we obtain 



log(l + Xr 2 ) = 



dw 2 

dZ2 



(3.68) 

(3.69) 
(3.70) 



Wi 



Z 2 



X 

V 



-I 



w 2 = log(l + Ar 2 ), z x 
2X(xx + yy) 



(a 2 — X(x 2 + y 2 + X(yx — xy) 2 ) 



dt 

y 2 

dt. 



(3.71) 



One can check that equation (3.71) transforms (3.6) to the form (2.3). 

Rewriting the first integrals I\ and I 2 in the integral form and identifying them 
in terms of the new variables, we have w\ = I\Z\ and w 2 = I 2 z 2 , which in turn also 
gives us a relationship between x and y with z\ and z 2 , respectively (after fixing 
the integration constants to be zero without loss of generality), that is, 



x = hzxy, 



1 + Xr A = e 



2 _ Jizi 



(3.72) 



Expressing x and y in terms of I\, I 2 , x and y (by utilizing the equation (3.68)) and 
substituting them in the expression for dz 2 , we obtain 



dz 2 



2 (Xa 2 r 2 - A(l + Xr 2 )(XI 2 + -^r 2 )^ 



/ 2 (1 + Ar 2 ) 

Now, from the expression for 1 + Ar 2 (vide equation (3.72)), we obtain 



dz 2 = — 



(a 2 



A 2 / 2 



h)e 



a 2 e- 21 ^ 



dt. 



Integrating the above equation we get 



1 



2Vh 



tan 



-2/ 2 + ( a 2 -A 2 / 2 +/ 2 )e 



-hz 2 



2Vh{{a 2 - X 2 I 2 + I 2 )e-^ - I 2 - a 2 e- 2I ^)h 



(3.73) 



(3.74) 



(3.75) 



where 73 is an integration constant which is nothing but the third integral of motion. 



In order to find the fourth integral, using dz\ 
dt to get 



dt 



dzi 



h 
2 



dz 2 



y 2 ((a 2 - A 2 / 2 + J 2 )e-^2 - I 2 - a 2 e -2/ 2Z2 )i 



and equation (3.74), we eliminate 



(3.76) 
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From equation (3.72), we get y 2 = (e l2Z2 
into equation (3.76) leads to 



(I 2 z 2 + 1) ~ 2 



l)/(X(I 2 zf + 1)), which on substitution 
(e /2Z2 - l)~ 1 dz 2 



_y 2 ((a 2 - X 2 I 2 + h)e~^ - I 2 - oP- e - 2I ^)h 
Now integrating (3.77) we get 

(I 2 - a 2 - A 2 / 2 ) + (a 2 - A 2 / 2 - I 2 )e^ 



tan 





1 , 


JlZl 


--tan- 



2A/i((a 2 - A 2 / 2 + I 2 )eJ^ - I 2 e 21 ^ - a 2 )i 



(3.77) 



(3.78) 



where I4 is the fourth integration constant. Now making use of these four integrals of 
motion, namely equations (3.68), (3.75) and (3.78), the general solution to equation 
(3.6) can be straightforwardly constructed. The resultant solution also agrees with 
the equation (5.40) of Chandrasckar et al. (2009), obtained through the modified 
PS approach, after a redefinition of integration constants. 

(vi) Example:9 Generalized linearizing transformation of Type-II 

To understand the generalized linearizing transformation let us start with the 
following system of coupled second order ODEs, 



k(xx — yy) + k 2 x , „ .. k(xy — yx) — k 2 y 
x+ v ; / y; ox + Xx = 0, y + 



(x 2 - y 2 ) 
The associated first integrals are 



(a; 2 - y 2 ) 



\y = 0. (3.79) 



h - 
h 



-XL 



k + \/^X(x + y) + x + y \ ^_ 
k — \/—X(x + y) + x + y 
k + V^A(x -y) + x-y Y 2v — 
k - V^A(a; — y) + x — y 



(3.80) 



Rewriting (3.80) in the form 
h 
h 



(k- 



-X(x + y) + x + y) dt 

fce -3v^At d 



1 



(k — s/—Xi{x — y) + x — y) dt 



'-X 
1 

7^ 



x + y 
k 

x-y 



+ 



k 



f-kt 



(3.81) 



and identifying the new variables, we get the linearizing transformation 

'■"-/ 
'•—I 



w i = H ; — ) e 



w 2 = (- 



k 

x-y 
k 



+ „) + *+*) 



dt, 



From the first integrals we have (after assuming the integration constants to be 
zero without loss of generality) 



wi = hzi, 



w 2 = I 2 z 2 . 



(3.83) 
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Utilizing (3.82) in (3.83), we get 



x{t) = - {hz l +I 2 z 2 )e-^ Xt 



^A 



y(t) = 2^ IlZl ~ l2Z ^ e ~ 



-A/ 



. (3.84) 



Substituting the expressions of x and y into equation (3.80), and solving the resul- 
tant equation for x and y, we get 



-A 



-{hz 1 +hz 2 )e V=Xt 



1 + e" 



-XI. 



1 _ e -V^A* 



'■(hzi - I 2 z 2 )e~ 



-A/ 



1 + e" 



- A/ 



1 



-A/ 



(3.85) 



Substituting the equations (3.84) and (3.85) into the expressions (3.82) for dz\ and 
dz 2l and integrating the resultant equation, we obtain 



Zl 



z 2 = 



1 f e 



h 

1 / ' < x 

i2 



-A 



+ (c 



-At 



-At 



/^(/s+tanh- 1 !- 
7 2 )(/ 4 +tanh" 1 [- 



-A< 



■]) 



-A/ 



■]) 



(3.86) 



where I3 and Z4 are the third and fourth integration constants, respectively. From 
(3.86) and (3.84), we can obtain the general solution for equation (3.79) straight- 
forwardly. 



4. Conclusion 

In this paper, we have studied the linearization of two coupled SNODEs. In partic- 
ular, we have introduced a new method of deriving linearizing transformations from 
the first integrals for the given equation. The procedure is simple and straightfor- 
ward. From our analysis we have demonstrated that one can have two wider classes 
of linearizing transformations, namely Class-A and Class-B, depending upon the 
nature of the independent variables. In Class-A category, the independent variables 
are the same, and we identified three types of linearizing transformations in which 
two of them arc new to the literature. On the other hand, in the Class-B cate- 
gory (the independent variables are different), we found six new types of linearizing 
transformations. We have explicitly demonstrated the method of deducing the lin- 
earizing transformations and the general solution for all of these cases with specific 
examples. However, in this paper we have restricted our attention on two aspects: 
(1) dependent variables are functions of (t, x, y) only, and (2) independent variables 
are not of the local form Zi = fi(t, x, y,x,y), i — 1, 2. Linearization under these two 
types require separate treatment and will be studied subsequently. The method 
proposed here can naturally be extended to any number of coupled second order 
ODEs and indeed one can derive a very wide class of linearizing transformations in 
these cases. 

The work of MS forms part of a research project sponsored by National Board 
for Higher Mathematics, Government of India. The work of ML forms part of a 
Department of Science and Technology, Government of India sponsored research 
project and is supported by a DST Ramanna Fellowship. 



Article submitted to Royal Society 



Integrability and Linearization 



19 



References 

Carinena, J. F., Ranada, M. F., Santander, M. & Senthilvelan, M. 2004 A non-linear os- 
cillator with quasi-harmonic behaviour: two- and n-dimensional oscillators Nonlinearity 
17, 1941-1963. 

Chandrasekar, V. K., Senthilvelan, M. & Lakshmanan, M. 2005 On the complete integra- 
bility and linearization of certain second order nonlinear ordinary differential equations. 
Proc. R. Soc. London A 461, 2451-2476. 

Chandrasekar, V. K., Senthilvelan, M. & Lakshmanan, M. 2006 On the complete inte- 
grability and linearization of nonlinear ordinary differential equations - Part II: Third 
order equations. Proc. R. Soc. London A 462, 1831-1852. 

Chandrasekar, V. K., Senthilvelan, M. & Lakshmanan, M. 2009 On the complete integra- 
bility and linearization of nonlinear ordinary differential equations - Part IV: coupled 
second order equations. Proc. R. Soc. London A 465, 609-629. 

Crampin, M., Martnez, E., & Sarlet, W. 1996 Linear connections for systems of second- 
order ordinary differential equations. Ann. Inst. H. Poincare, Phys. Theor. 65, 223-249. 

Duarte, L. G. S., Moreira, I. C. & Santos, F. C. 1994 Linearization under non-point 
transformations. J. Phys. A: Math. Gen. 27, L739-L743. 

Euler, N., Wolf, T., Leach, P. G. L. & Euler, M. 2003 Linearisable third-order ordinary 
differential equations and generalised Sundman transformations : The case X'" — 0. 
Acta. Appl. Math. 76, 89-115. 

Euler, N. & Euler, M. 2004 Sundman Symmetries of nonlinear second-order and third- 
order ordinary differential equations. J. Nonlinear Math. Phys. 11, 399-421. 

Fels, M. 1995 The equivalence problem for systems of second-order ordinary differential 
equations. Proc. London Math. Soc. 71, 221-240. 

Grossman, D.A. 2000 Torsion-free path geometries and integrable second order ode sys- 
tems. Selecta Math., New Ser. 6, 399442. 

Ibragimov, N. H. 1999 Elementary Lie Group Analysis and Ordinary Differential Equa- 
tions. New York: John Wiley & Sons. 

Merker, J. 2006 Characterization of the Newtonian Free Particle System in m > 2 Depen- 
dent Variables. Acta Appl. Math. 92, 125-207. 

Mahomed, F. M. & Leach, P. G. L. 1989 The Lie algebra sl(3, R) and linearization. 
Quaestiones Math. 12, 121-139. 

Mahomed, F. M. & Qadir, A. 2007 Linearization criteria for a system of second order 
quadratically semi-linear ordinary differential equations. Nonlinear Dynamics 48, 417- 
422. 

Qadir, A. 2007 Geometric linearization of ordinary differential equations. SIGMA 3, 103(1- 
7). 

Olver, P. J. 1995 Equivalence, invariants, and symmetry. Cambridge: Cambridge Univer- 
sity Press. 

Soh, C. W. & Mahomed, F. M. 2001 Linearization criteria for a system of second-order 
ordinary differential equations. Int. J. Non-Linear Mech. 36, 671-677. 

Sookmee, S. & Meleshko, S. 2008 Condition for linearization of a projectable system of 
two second order ordinary differential equations. J. Phys. A 41, 402001(1-7). 

Steeb, W. H. 1993 Invertible Point Transformations and Nonlinear Differential Equations. 
London: World Scientific. 



Article submitted to Royal Society 



